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The General Theory of Linear q-Difference Equations.* 

By R. D. Caemiohael. 



Introduction. 



By means of a transformation of the form z = (m-^x + rn^j^x -f- ^ 2 ), the 
system of functional equations 

Hi (£+1) = J>(^( z ) (* = 1, ■ ■ ■ ■ , «), 

in the n unknown functions S^z), . . . . , HJz), may be transformed into the 
system of difference equations 

n _ 

Glx + 1) = 2 VO #/*) (* = !,••••,«), 

or into the system of g-difference equations 

n _ 

(^faa;) = 2 V<b) G^jb) (• = 1, . . . . , n), 

according as the substitution z' = (az + o)J{cz + d) has one or two double points. 
To do this it is necessary to choose the transforming substitution so that in the 
first case the single double point is carried to infinity, and in the second case 
the two double points are carried to zero and to infinity respectively. 

The essential general properties of the solutions of linear difference equations 
are known.f 

The present paper is devoted to an investigation of the existence and prop- 
erties of solutions of linear ^-difference equations. 

In § 1, for the case when l^] =f= 1, I prove the existence of two fundamental 
systems of solutions, one of simple character at infinity and the other of simple 

* Read before the American Mathematical Society (Chicago), April 38, 1911. 

t See the- papers by Carmichael and by Birkhoff in Transactions of the American Mathematical Society, 
Vol. XII. References to the previous literature of difference equations will be found in these papers. 
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148 Oakmichael: General Theory of Linear q-Difference Equations. 

character at zero. These two systems of solutions are analogous to the two 
systems of solutions of difference equations whose existence I pointed out at the 
close of my difference-equation paper (loc. tit., pp. 133-134). 

In § 2 an investigation of the relations between these two fundamental 
systems of solutions leads to a theory analogous to the interesting Birkhoff 
characterization of the solutions of a system of difference equations (see §§ 5 
and 7 of the paper already referred to). 

In §3 I consider the exceptional case when |^|=1. A method is given 
for obtaining fundamental systems of solutions in explicit finite form. 

Rev. F. H. Jackson * has given a treatment of some questions connected 
with ^-difference equations. Reference should also be made to some of the 
general theorems on functional equations due to A. Gr6vy f and L. Leau, J from 
which some of the results of the present paper can be deduced. 

§ ]. Existence of two Fundamental Systems of Solutions when \q\^= !• 

ljet Glqx) = q«x" X A M {x) Gj(x) = qW X B^x) G } {x) (i = 1, . . . . , n) (1) 

j=i 1=1 

be a system of n first-order linear homogeneous ^-difference equations involving 
the n unknown functions G-^x), . . . . , G n (x) of the complex variable x, the 
known quantities entering into the equation being defined as follows: 

1. a, fi, q are constants and l^l =£ 1. 

2. The functions A v (x) and B^(x) are single-valued and 

A v {x) = A tj + A^x- 1 + 4jar+ + .... (i,j= 1, . . . ., n) \x\> B, 
B ij (x)=B v + B> v x + B<4x* + . . . . (»,/= 1, . . . ., n) \x\<r. 

3. The constants A tj and the constants B i} are such that the roots 

A lt . . ■ ., A n and B 1} , B n respectively of the characteristic equations at 

infinity and at zero, 

are finite and different from zero and verify the relations 

qiAi —ArfO, qTB t - BjjzO (t, j - 1, . . . . , n; i =£/), 

y being any positive or negative integer or zero. 

* Amebican Journal of Mathematics, Vol. XXXII, pp. 305-314. 
f Paris thesis, 1894. 
t Paris thesis, 1897. 
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When the roots of a characteristic equation satisfy the given relations, we 
shall say that they are of simple character. 
By means of the transformations 

Glx) = e^^^Glx) = e^°^^G' i {x) (i=l, ....,*), (2) 

where vj = log a/log q, system (1) may be replaced by two simpler systems, one 
of which has a = and the other of which has /3 = 0. Making these substitu- 
tions and reducing, we have 

n 

G t (qx) = 2 A ¥ (x)Gj(x) (i = 1, . . . • , »), (3) 

n 

G[{qx) = 2 B^x)G}(x) (i = 1, • • • • , n). (4) 

For the sake of simplicity in formulas it is convenient to reduce (3) and (4) 
to normal forms valid in the neighborhood of infinity and of zero respectively. 
We begin with system (3). 

Let US PUt Glx) = I aMx) (i = 1 , n), (5) 

5=1 

where oc^, i, /== 1, . . • ., n, are constants whose determinant |a#| is different 
from zero. Making this substitution in (3) and solving for Flqx), i = l, ....,«, 
one may write the result in the form 

n 

Flqx) = X A v {x) Fj(x) (»=1, »). 

)=\ 

where -the functions A tj (x), being linear combinations of the functions A^(x) with 
constant coefficients, may be expanded in powers of sc -1 for |sc| = i2. Prom the 
well-known theory of linear substitutions* and the fact that the roots of the 
characteristic equation at infinity are of simple character, it follows that a proper 
determination of the constants a^ will reduce the constant term in A v (x) to zero 
when i=f=j. The preceding system of equations may then be written in the 
form 

Flqx) = A, Fix) + I Kar 1 + o&ar" + . . . . )F,(x) (6) 

(* = 1 ,n) \x\>B. 

*Compare my diflerence-equation paper, loc. eit., pp. 136-127, where for a similar transformation the 
reduction is carried out in detail. 
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In order to remove the term in a; -1 from each of the non-diagonal coefficient 
functions of the second members, we may employ the transformation 

Fix) = Fix) + i &L Ffr), U = O, (i = 1, • • • • , n), (7) 

where the constants (3%, i=f=j, are to be determined. Clearly the determinant 
of the transformation is not identically zero. If we substitute the above values 
of F^x) in the preceding system of equations, and in the result write only those 
terms which do not involve x~ z , x~ s , , we have, when j x | = B, 

— n a _ 

Flqx) + 2&L Fj{qx) 

= AFlx) + AX^Flx) +|{(^ + • • • •)[*>) +1^^)]! 

where the accent in 2' denotes that in the summation the term for which j = * 
is to be omitted. 

In order to solve these equations for Flqx), i = 1, . . . . , n, in terms of 
Fix), i — 1, . . . ., n, we shall first find the cofactors of the determinant of the 
system, 



A = 



w qx 



= 1 + terms in x~ z , x~ 3 , | x \ = E. 



The cofactor Ay of the element in the i-th row and the^'-th column has evidently 
no constant term, except when i equals j, when there is a constant term 1. 
The term in a; -1 for i zfcj is clearly — fijijqx. For, if we cross out the i-th. row 
and the /-th column of A, there remain n — 2 elements with a constant term 1 ; 
and therefore the expansion of this minor can contain only one term in a; -1 , 
namely, that term which contains all the n — 2 factors 1 and whose remaining 
factor is not in the same rows or columns as these. Therefore this term is 
± fiji/qx ; and, since an even number of interchanges of rows and columns brings 
the elements @ v /qx and Px/qx to the place in the first row second column, and 
second row first column respectively, the negative sign must be chosen. In the 
cofactor A u there is no term in a; -1 ; for every term in the expansion containing 
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one factor cc -1 necessarily contains another. Hence the cofactors may be 
written 

** ~ ~~ qk + termS in X ~ 2 ' X ~ 3 ' "" ^ ^^ ' 



An = I + terms in a? -2 , as"" 8 , 



Solving the last system of equations by the aid of these cofactors, we 
readily obtain 

AFlqx) = (A, + **- + , • • • ) Fix) + £' \ ^ = Wv + ggg +,...[ I- (a .) 

(t = 1, ....,«). 

Since the roots of the characteristic equation at infinity are of simple character, 
it follows that it is always possible to determine /?^, i =fcj, so that (qA t — A f )fi v 
+ qa v = 0. If the values of (3 V so determined are substituted in the preceding 
system and each of the equations is then divided by A, the result takes the 
form 

(i=l,....,n) \x\lR. 
It is now easy to see that by means of the transformation 

Ffa) =/&)+% *£/&}, Y* = 0, (* = !,-• •-,»»), (9) 
system (8) goes over into the form 

ifoz) = (A + «« x- 1 + b u x-*)flx) + I (p v x-* + .... )flx) 

(i= 1, , «) | a; jiiB. 



j=i 



By repeated use of transformations similar to (7) and (9), the exponent of x in 
the coefficients of the transformation being increased by unity at each step, it is 
clear that our system of equations may be reduced to 

n 

flqx) = (A + Aix- 1 + . . . . + Afx-*)Mx) + jSJcfcr*- 1 + ..-. )/j(x) (10) 

(i=l,....,n) H>i2, 
where s is any positive integer. 
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By combining into one all the transformations thus far employed one 
readily proves that the change from, (1) to (10) may be effected at once by a trans- 
formation of the type 

G i ( B ) = e* to " WH ^^(i.)/X«) (.•=!,....,»), (11) 

where the it^ ( — J are polynomials in 1/x, the determinant of the transformation 

not being identically zero. 

In a similar way one may find polynomials P%{x) in x such that the trans- 
formation 

Glx) = e* loeq( ™2P«(xMx) (i=l, n) (12) 

is non-singular and carries (1) over into 

g i (qx) = (B i + B[x +....+ BFx°)glx) + I (^af +1 + . . . .)&(») (13) 

j—i 

(i= 1, , n) \x\lr, 

where a is any positive integer. It is hardly necessary to remark that the forms 
of the various preliminary transformations here are what those of the preceding 
case become when 1/a; is replaced by x. 

If A = q* f it is easy to determine constants c', c", . . . . , c (3s) such that 

f(x) = af(l + c'aT 1 + + c (s) aT s ) 

verifies the equation 

f{qx) = (A + 4'ar 1 + + 4 (s) a;- s )/(a;) + af(c (s+1) sr" 8 - 1 + + c (2s) x- 3s ). 

Evidently this may be written 

/ ,,(8+1) /»-»— 1 _|_ _1_ /i(2s) ™— 23 \ 

/(^) = (j. + .... + a^x-* + \ - ^ + ;;;; + ^ )/(*) = g ^)/H 

where /l(sc) has the constant term 1 and is analytic when | x | is greater than the 
greatest absolute value of a zero of 1 + dx~ x + .... + c (s) * _s - Consequently, if 
we write J^ = q\ there exist functions /l { (a;) and a constant R greater than or 
equal to B such that system (10) may be written in the form 

n 

flqx) - fMx)M*) = 2 M x )ffc) (» = 1, • • • • , «), (14) 

j'=i 

where each of the n 2 + n functions ^(as), ^(x) is analytic for | »):!$, the 
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expansion for 1^(05) beginning with a term in as -8-1 , and where the functions 
\{x) have each the constant term 1 and are such that for every i the equation 

fi 1) (qx)-q^ i (x)f< i \x) = 
has a solution in the form 

//"(a?) = UJx) = afi(l + C&- 1 + ....+ cfx~ s ). 

In a similar way it may be shown that if B % = q m < there exist functions 
l t (x) and a constant r different from zero and equal to or less than r such that 
system (13) may be written in the form 

n 

glqx) — q m aix)glx) = 2 ^x)g/x) (i = 1, . . . . , n), (15) 

3=1 

where each of the n* + n functions <p v (x), l t (x) is analytic for |a| = r, the ex- 
pansion for $ v (x) beginning with a term in x" +1 , and where the functions l t (x) 
have each the constant term 1 and are such that for every i the equation 

gf\qx)-q^l i (x)gf\x)-0 
has a solution in the form 

gf\x) = Vlx) = x m <{\ +d>x + .... + c^x"). 

Systems (14) and (15) are the normal forms of (1), the first being valid in the 
neighborhood of infinity and the second in the neighborhood of zero. To each of 
these normal forms the method of successive approximation is conveniently 
applicable. By means of this method we shall now prove the existence of a 
fundamental system of solutions of the equations in the normal form (14). 

Consider the set of systems each of n linear equations, all except the first 
being non-homogeneous : 

mqx)-q»<Ux)mx)=0, 

n 

ff\qx) - q^lx)ff\x) = 2 -*«, (x) /<»(*), 



7=1 



ff\qx) - q^{x)ff\x) = 2 *«(«)/FH 



1=1 



(16) 



We have already found a solution U t (x) of the first of these systems of equations ; 
we shall consider here the more general solution 

u i (x)=C i {x) U t (x) = C t (x) a*(l + ofari+ . . . . + c f) x -) (» = 1, . . . . , n ), (17) 
where O^x) is any function satisfying the functional equation G^qx) = C^x). 
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All the other systems have the general non-homogeneous form 

hlqx) — q^\{x)hlx) = ^(x) (i = 1, , n), (18) 

only a single unknown function h t (x) entering into the i-th. equation of the system 
t = 1, , n. Two formal solutions are readily obtained as follows : Let 

hi(x) =U i (x)h i (x) 

and substitute in the preceding equation. Dividing the result by 

U i (qx) = q*a i (x)U i (x), 
we have 

a system with the two formal solutions 

h (x) - — 2 ^fx) h(x) - 2 ^ ( J~"~ 1 ^ 

HX) ~ l«UWx)> hl{X) ~Zo U&T*) ' 

Hence two formal solutions of (18) are 

V.)=-r < (.)| | ^j. *M=VMl$g$. (w) 

Either of these may be used to obtain a sequence of formal solutions of (16). 
It will turn out that one of these sequences will lead to a (convergent) solution 
and the other to a divergent formal solution of (14). The first leads to a (con- 
vergent) solution if \q\ >• 1, the second if |gj<l. As the reasoning is entirely 
similar in the two cases, it is sufficient to carry it out for one. We shall suppose 
that \q\ > 1 and shall therefore use the first formal solution (19). 

Employing the notation 

and using for each system that particular formal solution which is obtained by 
adding the solution u^x) of the homogeneous equation to the solution derived 
from the first equation (19), we have 

ff\x) = Ut(x), 

n 

ff\x) = «,(») + SJ$$ v u,), 

n n n 

ff(x) = Ui{x) + S xi \ 2 hM + 2 4vStf(2 4fe«») } , 
j=i j=\ fc=i 
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By actual substitution one may readily verify that- a solution of system (14) is 

obtained in the limit functions (if they exist) of the sequences f^\ ff, , 

i = l, . . . . , n ; that is, in 

n n n 

fix) = u t (x) + SJ2, 4*yu>j) + Sm \ 2 4>y 5* (2 4>j k u k ) j 

3=1 j=l k=l 

n n n 

+ SJ& ^M 2 ^ jk S xk (Z ■*«!*,) I ] + . . . . (* = 1, . . . - , n). (20) 
j=i &=i i=i 

For a fixed j let us take O t (x) = 8^. We thus obtain a particular formal solution 
of (14), which we will denote by /#(#), i=l, . . . . , n. Making this substitution 
for every j, we obtain the n particular formal solutions 

M x )> M x )> ■■■■> M x ) (i = 1, ■•■-, n), 

where 

n 

fv(%) = hUA x ) + Sxi&ijU? + Sxi\^^ikS xk {^ M Uj)\ 

k=l 
n n 

+ S mt [24> ilt S xk \2 4, a SJ$ v U i )\] + .... (»,/ = I, ..-., n). (21) 
a=i i=i 

There exist constants R', M, c a and c 2 such that, for every i and/, 

\^{x)\<m\x\-^-\ Cl <:m(x)x-^\<:c„ \x\>&. 

Then the convergence of all the series in (21) and of the series (21) itself for 
|a;| = i2' will follow readily from similar convergence in the series 

^G&0+M££*-(A)}+ c> 

where 

We will suppose now that the arbitrary integer s has been chosen greater than 
twice the absolute value of any p. Considering the first term of (22), we have 



where q 1 denotes the sum of the last written series — a series which converges, 
since s>|/t^|. Prom this result it follows that the series afV /S^-^ty Z7J) is 
uniformly convergent for |sr |= R' and that its sum is less in absolute value than 
q 1 1 a; -8-1 1 . Moreover, each term of the series is analytic for | x | > R'; and, since 
20 



156 Cakmiohael: General Theory of Linear q-Difference Equations. 

the series converges uniformly throughout the region, the sum is analytic in the 
same region. 

Taking the second term of (22) and using the result of the last paragraph 
we have 



8 I § M 8 J Mc * M ' 8 \ Mn a x^-A = 1 Mnqi 






i[f77.C 2 



= gl i «fr*- 1 s o - 1 ^;; 2 ^,) i = «& I «^ g I » 



where g 2 denotes the sum of the last preceding series. From this we conclude 
that the third term in the series in (21) is in absolute value less than 
q^q % | art _2s-2 1 , and that it is the product of art by a function which is analytic 
for | a> | > 22'. 

Continuing this process we reach the result that every term in (21) is the 
product of art by a function which is analytic for \x\± R', and that the series is 
term by term less in absolute value than the series 

8<fy | art | + ft | art— 1 1 + q# t | tfr*-* | + qi q 2 q 8 \ art" 3 *- 3 | + . . . . , (23) 

where 

Mci _. _ ^ Mnc % 



qi ~^o c 1 |^ 8 + 1 -^|' qa ~i c 1 \4**~-*>\' a ~ 2 ' 3 ' 



Since q z , q 3 , . . . . is a decreasing sequence with the limit zero, it follows that 
series (23) converges when |a;| > 1. Hence there exists an 22", 22">1, R"> 22', 
such that the series obtained by multiplying (21) term by term by af~">- is uni- 
formly convergent and its terms are analytic for | x |=22". Hence all the functions 
x~^/ v (x), i,j = 1, ....,«, defined by (21) are analytic for |ar|> 22". 

Moreover, since the series in (21) is term by term less than the series (23), 
and since the limit as x approaches infinity of Uj(x)x~^ is 1, it follows that 

^=00 *-**'/«(*) = *» ft /=!,.• ■■,»)■ 

From this we see that the determinant of which the elements are f v (x) is not 
identically zero, and this is clearly a necessary and sufficient condition that 
the solutions f l3 {x) are linearly independent in the sense that no solution is 
expressible linearly in terms of the others, the multipliers being functions 
satisfying the equation 0(qx) = C{x). 
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If now we write 

GftXx) = e * lcsq(ri>+r,) 2 n ik (—) f v {x), n — log a/log q, (i, j = 1, . . . . , n), 
it is clear from (11) that we have n solutions 

&8\*), G&i*), ■•, <*${*) (j=l,....,n) 

of equation (1). Furthermore, since the transformation (11) is non-singular, it 
is evident that these solutions are linearly independent and constitute therefore 
a fundamental system of solutions of (1). On account of the properties of /#(«) 
it is clear that we may write 

anx)==e s^™ x », v $ Kx) iU = lt .... f n)i 

where Vffi{x) is analytic in the neighborhood of infinity. 

If an expression of the above form for G\J\x)-is substituted in (1), Vffi(x) 
being expanded as a power-series in 1/x, and if the constants are determined by 
direct reckoning, it will turn out that each of the n formal solutions obtained in 
this way is unique up to a constant multiplier. Consequently, if these constant 
multipliers are properly chosen, these formal solutions are the actual solutions 
whose existence has just been determined. 

This solution of (1) has been obtained by starting from the normal form (14) 
and employing the method of successive approximation. Starting from (15) and 
working in a similar way, it may be shown that we have also a fundamental 
system of solutions 

0#(*), <*${*), • • • • . &&{*) U = 1, n) 

which may be written in the form 

aif{x) = e> S9{M x^Vlf{x) ft/=l, ...., «), 

where 17°' (a;) is analytic in the neighborhood of infinity, and that this solution 
may be obtained by substituting in the equation as in the foregoing case and 
determining the constants in a direct way. The principal modification in the 
argument consists in using the second instead of the first formal solution (19) of 
the non-homogeneous equation (18). The discussion follows so closely along 
the lines of the preceding that it is unnecessary to repeat it. 

The argument above has been carried out on the assumption that | q \ >■ 1. 
If \q\ < 1, the only essential modification of the argument consists in inter- 
changing the roles of the two formal solutions (19) of (18). In this case the 
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second should be used for solutions about infinity and the first for solutions 
about zero. The results as to the nature of the solutions at these points still 
remain valid. 

The nature of each of the functions of the two fundamental systems of 
solutions obtained above is so far determined only in the neighborhood of infinity 
and of zero respectively. In order to determine their nature outside of these 
regions, we work as follows: If we consider the first two members of (1) and 

solve for G t (x) in terms of Giiqx), i = 1, , n, we may write the result 

in the form 

n 

0«OO = ?"" *~ a 2 A i} (qx) Q s (qx) (i = 1 ,...., n). (24) 
1=1 

Then from (1) itself and the preceding system, we have the relations 

n 

X 

k=l 
n 

2. 

k=l 



Gvfa) = «"*" 2 A ik (x) G kj (x) (»,/=!,...., n), (25) 

k=l 
n 

Gtt(x) = q~" x ~ a 2 A ik (qx) G kj (qx) (i, j = 1, . . . . , n), (26) 



where G i} (x) is to be identified with either Gffiipc) or G$(x). 

Let r lf r z , r 8 , . . . . be an infinite sequence of increasing numbers such that 

r i+i= l?K or r i+i= Iql^n, * = 1,2, , according as \q\ >1 or \q\ < 1. 

Consider the infinite system of circles of radii r lf r 2 , r 3 , . . . . and center at the 
origin, and assume that r x has been so chosen that the first two of these circles 
lie within the region in which V^f(x) is analytic. These circles divide the part 
of the plane external to the first one into an infinite system of circular rings, 
which we shall call fundamental regions. If we know the value of G$ (x) at 
every point in one of these rings, we can compute its value at any point in the 
next larger by means of (25) if \q\ > 1, and by means of (26) if |<z|<l- 
From this fact it is easy to determine the position of the singularities of G$(x). 
If fi=zq t v, where t is any positive or negative integer (not zero), we shall say 
that (i is externally or internally congruent to v according as \(i\ > \v\ or 

| fi | < | v | . Excluding the points zero and infinity, we see that the remaining 
singularities of each function G${x) are all included in the set of points exter- 
nally congruent to the singularities of A ik (x) or of A ik (x), i, k = 1, . . • • , n, 
according as | q \ >■ 1 or | q \ < 1 . 

In a similar way it is readily shown that the singularities of each function 

G ( ^ ) {x), other than zero and infinity, are all included in the set of points inter- 
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nally congruent to the singularities of A ik (x) or of A ik (x), i,k = l, , n, 

according as |<?|>1 or \q\ < 1. 

Furthermore, if A ik (x), i,Je = 1, . . . ., n, are rational functions, it is clear 
that the singularities (other than zero and infinity) of the functions in either 
solution consist entirely of poles. 

The principal results which we have obtained may be stated in the following 
theorem : 

Theorem I. The system of q-difference equations (1) or (24) has two funda- 
mental systems of solutions of the forms 

(?i->(x) = ei'°"" + V-7irW, « = jS§f, ft/=l,. ...,.). 

<?§>(*) = e^"- M ^7if(x) li,j= 1 „), 

where Yitf ) (x) and V$(x) are analytic in the neighborhood of infinity and of zero 
respectively. These solutions may be determined by substitution (a power-series being 
assumed for each of the functions V) and direct reckoning out of constants. Besides 
the points zero and infinity, all the singularities of Gj;f } (x) are included in the set of 
points internally congruent to the singularities of A ilc (x) or of A ik (x), i,kz=zl,. . . . ,n, 
according as \q\ ]> 1 or \q\ <C 1 ; and all the singularities of Glf(x) are included 
in the set of points externally congruent to the singularities of A ik (x) or of A ik (x), 
*, % = 1, . . . ., n, according as \q\ > 1 or \q\ <[ 1. The general solution G^x), 
% = l, . . . . , n, may be written in either of the forms 

n n 

G^x) = 2 OJT\x) G%\x) = 2 Ojf»{x) Glf{x) (i = l, . . . ., n), 

3=1 j=l 

where the functions G} a>) {x) and Gj 0) (x), j = 1, ...., n, satisfy the equation 
G(qx) = C(x) but are otherwise arbitrary. 

Further, if A ik (x), i,k = l,. . . . ,n, are rational functions of x, the singularities 
(other than zero and infinity) of the functions in each of the two fundamental solutions 
consist entirely of poles. 

Let us now consider a single ^-difference equation of the n-th order in the 
two forms 

E(q n x) + (q n xY A^x) H^- 1 x) + (q n x) a ((f' 1 x)« A z (as) S(q n ~ 2 x) 

+ ....+ (q n xY .... (q* x) a (qxY A n (as) E(x) = 0, ^ 

and 

H(q n x) + (q n xyB 1 (x)H(q n - 1 x) + {q n xy(q n - 1 xyB z (x)H(q n - z x) 

+ ....+ (q n xf .... tfxY(qxYB n (x) S(x) = 0, " (28) 
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where the unknown quantities entering into the equation are denned as follows: 

1. a, /?, q are constants and \q\ ^ 1. 

2. The functions A t (x) and B^x) are single-valued and 

A i (x) = A l + A' i x- 1 + Al'x- 2 + .... (t=l, ....,n) |a|>i2, 
£,(») = ^ + ^'8? + 5," «»+.... (» = 1, ....,») |a>|<r. 

3. The constants A % and the constants 5 4 are such that the characteristic 
equations at infinity and at zero, 

f + ^p"- 1 + Ap»- a +....+ 4_ x p + A. = 0, (29) 

p» + £ ip »-i + i? 2 p m - 2 + . . . . + ^-iP + B n = 0, (30) 

have their respective roots a x , . . . ., a n and b 1} .... 6 n of simple character. 
The transformations 

H{x) = 6 T * 5"(aj), 57 = log as/log <?, 

J( B ) = e frlfl " w+,fl 5>( B ), 
carry equations (27) and (28) respectively over into 

W(q n x) + A 1 (z)W(q»- 1 x) + ..-.+ 4,(05) IF (as) = 0, (31) 

H'(q n x) + B^x) H'{q n - X x) +....+ B n (x) H' (x) = 0. (32) 

It is evident that these two equations are equivalent respectively to the 
two systems 

E 1 (qx) = X\-Mx)\H j (x), E i (qx) = H i _ 1 (x) (i = 2, . . . ., n), (33) 

E[(qx) = l{-B j {x)\HJ(x), H>(qx) = Hl_ 1 (x) (• = 2, . . . ., n), (34) 

where H n (x) = H(x) and Hn(x) = H'(x). The characteristic equation of (33) 
at infinity takes the form (29), and the characteristic equation of (34) at zero 
takes the form (30). 

Equations (33) and (34) are of the form (3) and (4) respectively and conse- 
quently the preceding discussion is applicable to the former systems. Corre- 
sponding to the solved form (24) of (1) are analogous solved forms of (31) and 
(33) respectively; namely, 

if (as) + Mx) H{gx) + 2 2 (x) H{q*x) +....+ A n {x) H{q n x) = 0, (35) 

13 t (x) = H w (qx), (» = 1, • • - ., n— 1); H n (x) = 2 \— A n _ j:+1 (x)\ H^qx), (36) 

where the functions A i (x), i = l, » . . ., n, are readily expressed in terms of the 
functions A i (x), i = 1, , n. 
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If now we apply our previous results to the systems of equations and 
interpret them with reference to the single equation of the n-th order, we obtain 
the following theorem : 

Theorem II. The n-th order q-difference equation (27), (28) or (35) has two 
fundamental systems of solutions of the forms 

HP{x) = e> enM «™V$» (x) (• = !, ....,»), 

where Vi'" ) (x) and Vf\x) are analytic in the neighborhood of infinity and of zero 
respectively. These solutions may be determined by substitution (a power-series being 
assumed for each of the functions V) and direct reckoning out of constants. Besides 
the points zero and infinity, all the singularities of H la>) (x) are included in the set of 
points internally congruent to the singularities of A t (x) or of A i (x), i = 1, . . . . , n, 
according as \q\ > 1 or | q | < 1 ; and all the singularities of Hf 0) (x) are included 
in the set of points externally congruent to the singularities of A i (x) or of J. 4 (x), 

i = l, , n, according as \ q \ > 1 or \ q | < 1. The general solution H(x) may 

be written in either of the forms 

H(x) = £ Oj*(x) H^{x) = 1 Of (x) Hf>(x), 

9=1 j=l 

where the functions 0/ co) (x) and Gj 0) (x), j = 1, , n, satisfy the equation 

C(qx) = C(x) but are otherwise arbitrary. 

Ftirther, if A i (x), i = 1, . . . ., n, are rational functions of x, the singularities 
(other than zero and infinity) of the functions in each of the two fundamental solutions 
consist entirely of poles. 

Remaek. In view of the preceding work it is obvious that the results of 
this section may be obtained by direct reckoning out of the solutions, by the 
power-series method suggested above, together with appropriate convergence 
proofs. Such a method, however, is not available in the case of difference 
equations. Consequently it has seemed best to employ the method of successive 
approximation (already used by me for difference equations), so as to develop 
the two parts of this general theory along the same lines and by the same 
methods. The matrix method used by Birkhoff for difference equations is also 
applicable to ^-difference equations ; it furnishes therefore a second single method 
by which both cases may be treated. 
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§ 2. Relations between the two Fundamental Solutions. 

The question of the relations existing between the two fundamental solutions 
of system (1) furnishes an interesting problem. It is clear that either of these 
solutions may be expressed in terms of the other. Thus we have 



Gff(x) = X G ik {x) GQ\x) (ij = 1 ,n), (37) 

Qip{x) = 1 C ik (x) G$(x) (i, J = l,....,n), (38) 



fc=i 



where G ik (x) and G ik (x) are functions satisfying the equation C(qx) = C(x). 
Our problem reduces, then, to that of the determination of the nature of these 
functions G tk (x) and G ik (x). In the general case they are of a complicated 
character. But in the most interesting and important case, namely, that in which 
the coefficients of the original equation are polynomials (or even when they are 
merely rational functions), the result takes a simple and elegant form. 

We suppose that system (1), subject to the previously assigned conditions, 
is also such that it may be written in the form 

n 

Glqx) = 2 P ik {x) a k (x) (* = 1, . . . . , n), (39) 

where the coefficients P ik (x) are polynomials in x of degree a or less. The 
constant (3 which enters into system (1) is in this case zero. 

According as | q | > 1 or | q \ < 1, the functions G$(x) or Glf\x) are 
analytic except at zero and infinity. We will consider first the case in which 
\q\ > 1. If (37) is solved for G ik (x), i, h = 1, . . . ., n, in terms of Gff(x) and 
a i} (x), it is readily seen that every function C ik (x) is analytic except perhaps 
at zero, innnity, the singularities of G^fa), and the zeros of the function 
defined by the* determinant | G^p (x) | . But on account of the properties Of the 
functions G^^x) and of the form of the determinant (see Theorem I), it is clear 
that there exists a fundamental region (sufficiently remote from the origin) in 

which the functions G$\x), i,j=zl, , n, have no singularity and the 

function | Gffi (x) | has no zero. Hence the functions G ik (x), i,h = l, . . . ., n, 
are analytic in this fundamental region; and therefore they are analytic through- 
out the plane, except at zero and infinity. 

In case \q\ <^ 1, we may argue in a similar way to show that the functions 
G ik {x), i,Jc = l, . • • • , n, are analytic throughout the plane, except at zero and 
infinity. 
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It is clear that if the coefficients P ik (x) in (39) are restricted to be rational 
(but not necessarily polynomial), a discussion analogous to that in the foregoing 
case can be readily carried out. If \q\ >1, the functions G ik (x), away from 
zero and infinity, will be analytic except for poles; and if \q\ < 1, the functions 
C ik (x) will have this property. 

These results may be stated in the following theorem : 

Theorem III. Let a system of q-difference equations (1) be written in the form 

G t (qx) = 1 P tk (x) G k {x) {% = 1, . . . . , n) ; 

and let the two simple fundamental solutions be connected by the relations 

G$(x) = I C (k (x)G^ (x), C ik (qx) = G ik {x) (i, j, k = 1, . . . . , n), 

where r, s = 0, oo or qo , according as | q | > 1 or \ q \ < 1. 7%e» : 

// <Ae coefficients P ik (x) are polynomials in x, the functions G ik (x) are analytic 
except at zero and infinity ; 

If the coefficients P ik (x) are rational in x, the functions G ik (x), away from zero 
and infinity, are analytic except for poles. 

An inverse to some of the preceding results may be stated in the form of 
the following theorem : 

Theorem IV. Let Gf}(x) and GfiXx), i,j=l, , n, be two sets of n z 

single-valued functions which, away from zero and infinity, are analytic except for 
poles, and which have the property 

J= 0ffO* ^ 10 " 9<,5+,) *-™ = V< $> I F $ I * ° i m <- ™J * int eger, i * /, 

J=. oo ( i?\^° Sqi ^ ) *-*' = V\l\ I y\T 1*0; in- n* integer, i ±J, 

where a and /? are integers (positive, negative or zero), ri = log a;/ log q, the two sets 
of functions being connected by the relations 

G\J(x) = I G ik (x) Gff(x), C ik (qx) = G ik (x), (i,j,k=l,...., n). 

Then the lico sets of functions Gff(x) and G['j ) (x) are the two simple fundamental 
solutions of a system of q-difference equations 

G t (qx) = I P ik (x) G k (x) (i=l ,n), 

in which the coefficients P ik (x) are rational in x. 
21 
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The proof of the theorem requires the determination of the existence and 
the properties of functions P ik (x) such that 



G$ (qx) = 2 P ik (x) Gftx) (i, j = 1, . . • • , «), (40) 
G$\qx) = I,P Vc (x)G<g\x) (»,/=!, ....,n). (41) 



4-1 



Each of the two systems of equations may be solved for Py^x). The two repre- 
sentations of P ik (x) must be consistent, in view of the relations connecting 
G$(x) and G§\x). From either of them it is clear that the functions P tk (x), 
away from zero and infinity, are single-valued and analytic except for poles. 

If (40) is solved for P tk (x) in terms of G§\x) and Gf](qx), i,j,k=l,....,n, 
and the limit for x = is taken, then, in view of the first limit in the 
hypothesis of the theorem, we have 

lim p , x _„_ |F«*>| 
x = ^ W x 1 7(0)i" i 

where | F<*' | is the determinant with T™ in the r-th row and s-th column and 
\V? s k) \ is what \V { ®\ becomes when the i-th row is replaced by q m iV§\ 

j = 1, , n. From this it follows that P ik (x) has at x = a zero of order /? 

(or a pole of order — /? if /? is negative). Furthermore, the ^-difference 
equation has at zero the characteristic equation 

I yak) i 

I Pik ^*P | = 0, P ik = jy* 



(0)1 



A direct reckoning will show that the roots of this equation are p = q m s, 
3 — 1, , n. 

In like manner we may employ (41) and the second limit in the hypothesis 
of the theorem to show that P ik (x) has at infinity a pole of order a (or a zero of 
order — a if a is negative) and that the roots of the characteristic equation at 
infinity are p = q*>, j = 1, , n. 

From the fact that P ik (x) is analytic except for poles we conclude that it is 
rational in x. From the characteristic equations at zero and infinity respectively, 
we conclude that the ^-difference equation of the theorem has solutions of the 
requisite property at each of these points. This completes the demonstration 
of the theorem as stated. 
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§3. The Case when \q\—\. 

When \q\ = 1, the theory of the ^--difference equation is essentially different 
from that when |gi=£l, previously considered. In treating the former we 
distinguish two cases, according as there is or is not a positive integer n such 
that q n = 1. 

Case I. We will suppose first that there exists no integer n such that q n =l. 
That there is in this case in general no analytic solution (not identically zero) 
is easily seen from the consideration of a first-order equation with constant 
coefficients : 

f(qx) = af(x), o>l. 

We have evidently 

f{fx) = a/iq^x) = a?f(q°-*x) = .... = a a f(x). 

For any x it is clear that there exists an a greater than any preassigned M 
and such that q*x is less than any preassigned distance e from x. When f(x) 
is not equal to zero at the given point a; the quotient of the functional values at 
q*x and x, that is f(q a x)Jf(x), can therefore be made as large as one pleases, 
although q a x and x are separated by a less distance than e. Hence the equation 
f(qx) = af (x) has no continuous solution except /(a;) = 0. 

Case II. We will suppose now that q is a root of unity and that n is the 
least integer such that q n = 1. In the present case it is more convenient to 
consider from the outset a single equation than a system of equations, for 
evidently any single equation of any order can be written in the form 

*i(«0/GpO + 4*(*)/(fl"«) +••••+ *n(x)f(q n x) = o. (42) 

That solutions (not identically zero) do not exist when the 4-'s are unrestricted 
is readily shown ; for we have the following system of equations : 

M»)ffa) + M*)f(tf*) +•■••+ *«(*)/(2"») = 0, 

4*(g»)/(«a») + Mq*)ftf*) +•■••+ ^n-i(^) /(«"») = o, 
*-i(s"«0/(2*) + Mq*x)f{q 2 *) +••••+ +-,(2 8 «) fdf*) - o, - (43) 
•• 
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This system is consistent only if its determinant is zero; that is, if 
4-1 (x) 4>z(x) '"Ma) 



^n (?») 4l (?») 4>n-l iSPO) 

4«-i (? 2 as) -kn itfx) .... 4-„- 2 (q 2 «) 
^ a C?"- 1 *) V (q™- 1 as) .... ^ (j"- 1 a) 



= 0. (44) 



We therefore suppose that the coefficients 1^(2) are connected by the relation (44). 

So far the only restriction on q is that it is a primitive n-th root of unity. 
The numbers q, q % , q 8 , . . . ., q n ~ l , q n are all the n-th roots of unity, and hence 
this set of powers contains the primitive n-th root of unity of smallest argument. 
Without loss of generality we may suppose that this root is q itself; and this 
we shall do. Then the points q, q 2 , . . . . , q n are situated on the unit circle at 
equal intervals, in such wise that in making a positive circuit along the circum- 
ference the points are encountered in the order written. If rays are drawn from 
the origin through these n points, the plane is divided into n equal sectors ; and 
these we shall call a set of fundamental regions of the plane. If these rays are 
all rotated about the origin through the same angle, then in their new position 
they will define another set of fundamental regions. 

We have seen that a necessary condition for the existence of a solution 
(not identically zero) is that the determinant in (44) shall be of rank n — 1 at 
most. We suppose in general that this determinant is of rank n — a, where a 
is a number of the sequence 1, 2, , n — 1. Then, among the quantities 

f{qx),f(q*x), ...., /(a**), 

there is a set a in number in terms of which every other can be expressed 
linearly. Hence there are a fundamental regions in which f(x) may be assigned 
at will ; it is then determined throughout the entire plane. On account of this 
fact we shall say that the order of the equation is a. 

When the equation is of order 1, in accordance with this definition, there 
is a single relation among the coeflBcients; otherwise there is more than one. 
The general case is therefore that in which the order is 1. Since in this case 
a solution is determined when its value is assigned in a single fundamental 
region, it is clear that the most general solution is of the form G(x)f 1 (x), 
where /i(a;) is any particular solution (not identically zero) and G(x) is an 
arbitrary function satisfying the equation C(qx) = C{x). 
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In the case of a first-order equation it is easy to obtain a particular solution 
in explicit finite form. The ordinary theory of linear algebraic equations applied 
to the system (43) yields us the values of the ratios 

f(qx):f(q*x):....:f(q»x); 

and hence we have an equation of the form 

f{qx) = +(x)f(x), (45) 

where f(x) is the same as the f(x) in (42). To find a solution of this equation, 
we proceed thus : 
The product 

f(q n x)f(q»-*x)....f(qx) 

is unchanged by the substitution of qx for x; and hence its value is a function 
satisfying the relation C(qx) = C(x). Consequently there exists a particular 
solution /j (x) for which 

/i(fl**)/i(S*- 1 »)-----/x(?*) = l. 
But 

Mq i+1 x) = ^(q'x)f 1 (q i x). 

Making repeated substitutions from this formula into the preceding equation, 
we have 

ibgr^tyisr'x)]'.... ["Kg*)]*- 1 [*(»)]» [/i (»)]" = l. 
Hence 

/i0») = {^{q n - 1 x)^{q n - i x)Y....\^{qx)Y- 1 i^{x)T\- Vn 

is a particular solution of (45) and hence of (43); the general solution is 
C(x)f 1 (x), where C(x) satisfies the equation G(qx) = C(x) but is otherwise 
arbitrary. 

The above discussion suggests a general method for solving (42) when the 
rank of its determinant in (44) is n — a. In a form most convenient for this case 
the method may be described as follows: To any n — a independent equations 
of the set (43) adjoin other equations subject to the following conditions : 
1) The n equations of the resulting set are independent and consistent. 2) The 
left members of the a new equations are invariant under the substitution of qx 

for x. They are linear in f(q n x), ,/(»)• 3) The second members of these 

a equations are arbitrary functions of x satisfying the relation C(qx) = G(x). 
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One of the possible forms for these equations is 

Mq n - 1 x)f(q n - 1 x) + . . . . +X(spe)f(qz) + *(*)/(») = <?(*), 

where % (x) is any conveniently chosen function of x. 

If the system is solved algebraically for f(q n x), , /(«), the resulting 

expression for f(x) will contain linearly a arbitrary functions which satisfy the 
equation G(qx) = G (x). And hence this gives the general solution of (42), 
since it is of order a. 

In carrying out this process it is important to secure that the arbitrary 
functions shall enter into the solution independently ; that is, in such way that 
they can not be combined in the solution into fewer functions of the same type. 

We shall illustrate the above method by a treatment of the special case 
in which q = a = cube root of unity of least argument and a = 2. It is easy 
to show that the general equation for this case may be written in the form 

f(a z x) + ^(x)^(ax)f(ax) -\-^(x)f(x) = 0, ^ (x)^(qx)^(o % x) = 1. 

In accordance with our general method, if ^ (x) =fc 1,* we may adjoin the 

equations 

f(tfx)+f(c>x)+f{x) = C 1 (x), 

a 2 xf(o 2 x) + axf(ax) +/(*) = G s (x). 
The value of f(x) is obtained algebraically from the last three equations. 

Indiana Umivbbsitt, July, 1911. 

* If tyix) = 1, a fundamental system of solutions is 1 and x. 



